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COORDINATE  TRANSFORMATION  FOR  PHASED  ARRAY  ANTENNA 
BEAM  STEERING  USING  GPS  AND  SHIP’S  MOTION  DATA 


Background 

In  order  to  steer  the  pencil  beam  of  a  phased  array  antenna  to  a  target,  the  target’s 
position  needs  to  be  known.  This  info  can  come  from  a  tracking  radar  or  from  GPS  data 
of  both  the  target  and  antenna  position.  When  the  antenna  is  installed  on  a  moving 
platform  such  as  a  ship,  then  this  antenna  is  sometimes  mounted  on  a  stabilized  platform. 
Such  a  stabilized  platform  is  not  required  if  the  antenna  can  be  stabilized  electronically. 

In  the  following,  the  algorithms  required  to  steer  the  pencil  beam  of  a  phased 
array  antenna  towards  a  moving  target  are  given.  The  antenna  is  mounted  on  a  ship.  Input 
are  GPS  data  of  the  position  of  the  antenna,  the  GPS  data  of  the  target,  and  roll,  pitch  and 
bearing  data  of  the  ship.  The  programs  given  were  written  in  MATLAB  and  can  be  run 
on  any  PC.. 

Conversion  of  the  geocentric  GPS  data  to  longitude,  latitude  and  height 

The  GPS  data  can  be  given  in  terms  of  longitude,  latitude,  and  height  above  sea 
level,  or  in  Geocentric  Coordinates.  If  the  data  are  given  in  Geocentric,  then  these  can  be 
transformed  to  longitude,  latitude  and  height  above  sea  level  as  follows; 

Figure  1  shows  the  convention  for  the  geocentric  convention. 

The  coordinate  system  is  of  the  Right  Hand  Cartesian  type  with  the  apex  in  the  center  of 
the  earth,  the  x-axis  passing  thru  the  intersect  point  of  the  equator  and  the  Greenwich 
meridian,  and  the  z-axis  passing  thru  the  north  pole. 

Longitude  LO,  latitude  LA,  and  the  height  H  above  the  sea  surface  at  the  equator  are 
calculated  from  the  Geocentric  x  y  z  coordinates  ; 

sin(LA)  =  z  /  r  ,  with  r  =  sqrt(  x^  +  y^  +  z^) 
tg(LO)  =  y/x 


H  =  r- 4E7  /  (2*71:)  ,  with  r  and  H  in  meters  . 

This  transformation  is  programmed  in  the  Matlab  program  gztolola.m  (see  appendix), 
with  x  y  z  the  inputs  and  LA,  LO  and  H  the  outputs. 

Rotation  of  the  coordinate  system 

The  coordinate  system  with  unity  vectors  el  ,  e2  ,  e3  in  the  x,y,  and  z  direction 
respectively,  is  rotated  resulting  in  a  new  coordinate  system  with  unity  vectors  epl  ,  ep2, 
and  ep3  . 

Then,  [1] 

epl  =all*el  +  al2*e2  +  al3*e3 
ep2  =  a21*el  +  a22*e2  +  a23*e3 
ep3=a31*el  +  a32*e2  +  a33*e3 
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where 


all  =cos(epl,el)  al2  =  cos(epl,e2)  al3  =  cos(epl,e3) 
al2  =  cos(ep2,el)  a22  =  cos(ep2,e2)  a23  =cos(ep2,e3) 

al3  =  cos(ep3,el)  a23  =  cos(ep3,e2)  a33  =  cos(ep3,e3) 

The  xl  x2  x3  coordinates  of  the  system  el  e2  e3  are  transformed  into  the  system  epl  ep2 
ep3  to  yield  the  new  coordinates  xpl  xp2  xp3  in  the  rotated  system: 

xpl  =  all*xl  +  al2*x2  +  al3*x3 

xp2  =  a21*xl  +  a22*x2  +  a23*x3 

xp3  =  a31’'‘xl  +  a32*x2  +  a33*x3 


Rotation  around  xl  -  axis: 

As  is  illustrated  in  Figure  2,  the  xl  x2  x3  coordinate  system  is  rotated  counterclockwise 
when  visioning  along  the  xl-axis  in  the  negative  direction ,  counting  the  angle  a  positive 
when  rotating  in  the  counter  clockwise  direction.  Then,  as  can  be  read  from  Figure  2, 

al  l=cos(epl,el)=l  ;  al2=cos(epl,e2)=0;  al3=cos(epl,e3)=0; 
a2  l=cos(ep2,el)=0;  a22=cos(ep2,e2)=cos(a);  a23=cos(ep2,e3)=cos(90-a)=sin(a) 

a3  l=cos(ep3,el)=0;  a32=cos(ep3,e2)=cos(90+a)  =  -sin(a);  a33=cos(ep3,e3)=cos(a) 

This  relationship  can  also  be  expressed  in  the  form  of  Innerproducts,  remembering  that 
the  Innerproduct  of  two  unity  vectors  equals  the  cosine  of  the  angle  between  them. 
Therefore,  al  1  =  cos(angle  between  epl  and  el)  =  inner  product  of  epl  and  el.  The 
inner  product  of  two  vectors  x=[xl  x2  x3]  and  y=[yl  y2  y3]  is  xl*yl+x2*y2+x3*y3  . 
We  may  write  this  as  “inrprod(x,y)  =  xl*yl+x2*y2+x3*y3  “  . 

However,  epl  and  el  must  be  in  the  same  coordinate  system,  epl  appears  in  the  original 
coordinate  system  as  eppl,  where  (see  Figure  2) 

eppl  =[1  0  0]  epp2  =  [0  cos(a)  sin(a)]  epp3=[0  -sin(a)  cos(a)] 
and 

el=[l  0  0]  e2  =  [0  1  0]  e3=[0  0  1] 
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then , 

all  =  cos(eppl,el)  =  inrprod(eppl,el)  =  1*1  +0*0  +0*0=1 
al2  =cos(eppl,e2)  =  inq3rod(eppl,e2)=  1*0  +0*1  +0*0  =0 
al3  =  cos(eppl,e3)  =  1*0  +  0*0  +0*1  =0 
a21  =  cos(epp2,el)  =  0*1  +  cos(a)  *0  +  sin(a)  *0  =0 
a22  =  cos(epp2,e2)  =0*0  +cos(a)*l  +  sin(a)  *0  =cos(a) 


Now,  with  all  the  anm  given,  the  new  coordinates  xpl  xp2  xp2  can  be  calculated  from  xl 
x2  x3  . 

The  Matlab  function  xp=rlx(x,a)  yields  xp,  where  x  is  the  vector  [xl  x2  x3]  in 
the  original  coordinate  system  and  xp  is  the  vector  [xpl  xp2  xp3]  in  the  coordinate 
system  after  rotation  counterclockwise  around  the  xl  axis  ,  counting  the  angle  positive 
looking  along  the  positive  xl-axis  towards  the  origin. 

Rotation  around  x2  -  axis 

Similarly,  for  rotation  around  the  x2  -  axis  the  vector  epp  is  found  (Figurel  2A): 

eppl=[cos(a)  0  -sin(a)]  epp2=[0  1  0]  epp3=[sin(a)  0  cos(a)] 

The  input  vector  x  is  transformed  into  the  vector  xp  by  rotation  around  the  x2 
axis  by  angle  a  by  the  MATLAB  function  xp  =  rly(x,a),  see  appendix. 

Rotation  around  x3  -  axis 

For  rotation  around  x3-axis,  the  vector  epp  is  (Figure  2B): 

eppl=[cos(a)  sin(a)  0]  epp2=[-sin(a)  cos(a)  0]  epp3=[0  0  1] 

The  vector  X  is  transformed  to  the  vector  xp  by  rotation  around  the  x3  axis  by  the 
MATLAB  function  xp  =  rlz(x,a) ,  see  appendix. 


Example 

Assume  that  a  vector  x  is  given  in  a  coordinate  system  xx,yy,zz  with  the  xx-axis 
directed  towards  the  east ,  the  yy-axis  towards  the  north  and  the  zz-axis  up.  To  be 
calculated  is  the  vector  xp  referenced  to  a  coordinate  system  centered  on  a  ship,  with  the 
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x-axis  parallel  to  the  length  axis  of  the  ship,  positive  towards  the  bow,  the  y-axis  parallel 
to  the  deck,  positive  towards  port,  and  the  z-axis  positive  up. 

The  coordinate  system  xx  yy  in  which  the  vector  x  is  given  is  rotated  around  the 
zz  axis  and  then  around  the  yy  axis  until  the  axis  of  the  rotated  system  falls  on  the 
respective  axis  of  the  ship’s  coordinate  system. 

The  transformation  is  performed  in  two  steps.  First,  the  vector  x  is  transformed 
into  a  coordinate  system  rotated  by  the  angle  TtH  -  beara  around  the  z-axis  (see  Figure 
3).  Now  the  projection  of  the  ship’s  x-axis  onto  the  xx  yy  plane  lines  up  with  the  x-axis 
of  the  rotated  system.  This  is  performed  by 

xp  =  rlz(x,7r/2-beara) 

The  second  step  is  to  rotate  the  rotated  system  around  its  y  -  axis  by  the  angle 
pitcha,  thereby  lining  up  the  x-  axis  of  this  rotated  system  with  the  ship’s  x-axis.  Now  all 
of  the  axis  of  the  ship’s  system  are  lined  up  with  the  axis  of  the  system  after  the  second 
rotation,  and  therefore,  the  resulting  vector  xp  is  in  reference  to  the  ship’s  system.  The 
commands  are; 


x=xp 

xp  =  rly(x, pitcha) 

Be  it  is  assumed  that  the  ship’s  bearing  is  45  degrees,  that  is,  alf  =  90  -  45  degrrees.  The 
pitch  of  the  ship  be  45  degrees  also.  Then  a  vector  x=[l  1  sqrt(2)] ,  when  transformed 
to  the  ship’s  coordinate  system,  should  be  in  the  length  axis  of  the  ship  and  have  a  length 
of  2,  i.e.  we  should  find  that  xp  =  [2  0  0], 

First  transformation:  roll  around  z  =  axis: 
x=[l  1  sqrt(2)]; 
xp  =  rlz(x,7r/4); 

Second  transformation:  roll  around  y  -  axis  to  compensate  pitch: 
x  =  xp; 

xp=  rly(x,  -ttM) 
results  in 

xp=[2  0  0] ,  which  is  correct. 

Note  that  the  negative  of  the  pitch  angle  had  to  be  used,  because  rolling  around  the 
y-axis,  looking  towards  the  origin,  the  original  coordinate  system  has  to  be  turned 
clockwise  towards  the  pitch  angle  of  the  ship.  However,  the  convention  was  that  the 
angle  is  counted  positive  for  counter  clockwise  (ccw)  rotation;  hence,  the  angle  is 
negative  for  clockwise  (cw)  rotation. 
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Transformation  of  the  Roll  angle  of  the  ship 

The  roll  angle  “rolla”  of  the  ship  is  defined  as  positive,  if  the  roll  around  the  length  axis  of 
the  ship  as  seen  in  direction  to  the  bow  is  cw.  Therefore,  with  this  convention,  when  seen 
towards  the  origin,  the  rotation  is  ccw.  The  transformation  is  then 


x=  xp 

xp  =  rlx(x,  rolla) 

Transformation  due  to  antenna  position 

So  far,  the  antenna  was  assumed  to  be  located  such  that  its  boresite  was  parallel  to 
the  ship’s  length  axis,  radiating  forward.  Let  “antposaz”  be  the  angle  by  which  the 
antenna  is  rotated  around  an  axis  which  is  perpendicular  to  the  ship’s  deck,  antposaz 
being  positive  for  cw  rotation  as  seen  from  on  top.  Since  the  rotation  is  cw,  the  angle 
antposaz  must  be  entered  negative  into  rlz.  Therefore,  the  commands  to  take  care  of  the 
position  of  the  antenna  are 


x  =  xp; 

xp  =  rlz(x, -antposaz) 

Transformation  due  to  tilt  of  the  antenna  boresite 

The  tilt  angle  “antposel”  of  the  antenna  boresite  is  defined  to  be  positive  when  tilted 
upwards.  The  coordinates  in  reference  to  the  antenna  are  shown  in  Figure  5.  Upward  tilt 
means  cw  rotation  around  the  y-axis  ,  looking  towards  the  origin.  Since  the  roll  algorithm 
is  defined  as  positive  for  ccw  rotation,  the  angle  antposel  must  be  entered  negative  into 
the  function  rly.  Therefore, 


x=xp; 

xp  =  rly(x, -antposel) 

Now  we  have  the  original  vector  x  given  in  coordinates  referenced  to  the  antenna,  with 
the  center  of  the  coordinate  system  in  the  center  of  the  antenna  and  with  the  x-axis  lined 
up  with  the  boresite  and  pointing  in  the  direction  of  the  radiation.  The  y-axis  is  pointing 
towards  the  left,  as  seen  from  on  top,  looking  down  on  the  antenna. 

With  xp  =  [xp(l)  xp(2)  xp(3)], 

The  azimuth  steering  angle  “azim”  of  the  beam  is  then 

azim  =  atan2(xp(2),xp(l)) 

It  is  positive  for  steering  the  beam  to  the  left,  as  seen  in  direction  of  propagation. 

The  elevation  angle  “elev”  of  the  beam  is 

elev=  atan2(xp(3),sqrt(xp(l)^  +  xp(2)^)) 

It  is  positive  for  the  beam  pointing  up. 
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Example; 


Latitude 

Longitude 

Heigth 

(degrees) 

(degrees) 

(m) 

Target; 

gpsplati  =  29 

gpsplong  =  -85.3 

gpsph=150 

Ship 

gpsslati=29 

gpsslong=-85.5 

gpssh=5 

Bearing  of  ship;  beara=  190  degrees 
Roll  anngle  of  ship  rolla=-5  degrees 
Pitchangle  of  ship  pitcha=5  degrees 

Azimuth  position  of  antenna  on  ship  antposaz=270  degrees  (port  side) 

Tilt  of  antenna  boresite  is  15  degrees  up;  antposel=15  degrees 

Figure  6  shows  the  geometry.  Both  target  and  ship  are  at  the  same  latitude;  the 
target  is  east  of  the  ship  by  0.2  degrees,  i  .e.,  12  minutes  of  arc,  or  12  nautical  miles,  or 
22.2  km.  At  that  distance  and  the  target  height  of  150m,  the  target  is  practically  on  the 
horizon,  as  seen  from  the  ship.  The  data  above  were  entered  into  testd.m  (see  appendix). 
The  commands 
testd 
laurant 

result  in 

azim=10.0  degrees 
elev=-8.6  degrees 

for  the  beam  steering  angles  of  the  antenna. 
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Summary  of  symbols  and  conventions 

X  =  The  vector  pointing  from  the  ship  antenna  to  the  target.  Its  length  is  equal  to  the 
distance  from  the  antenna  to  the  target.  Its  coordinate  system  is  Cartesian,  with 
the  x-axis  pointing  to  the  east,  the  y-axis  to  the  north  and  the  z-axis  up.  This 
vector  is 

transformed  by  several  rotational  transformations  to  the  coordinate  system 
centered  at  the  antenna. 

beara  =  the  ship’s  bearing  angle,  positive  going  cw  starting  at  true  north 
pitcha  =  the  ship’s  pitch  angle,  positive  with  the  bow  up 

rolla  =  the  ship’s  roll  angle,  positive  rolling  cw  as  seen  looking  towards  the  bow 

antposaz  =  azimuth  position  of  the  antenna  on  the  ship,  cw  starting  from  the  bow. 

antposaz=0  for  the  boresite  parallel  to  the  length  axis,  radiation  in 
forward  direction 

antposel  =  elevation  angle  of  boresite,  referenced  to  the  deck  of  the  ship,  positive  up. 


azim  =  azimuth  pointing  angle  of  beam,  offset  from  boresite  in  azimuth.  Pointing  to  the 
left,  as  seen  indirection  of  propagation  renders  positive  angles. 


elev  =  elevation  pointing  angle  of  beam,  positive  up. 


[1]  Hutte,  Des  Ingenieurs  Taschenbuch ,  28  th  edition,  p.  1 15 
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Rotation  Around  X3  Axis 
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Geometry  of  Example 


APPENDIX 


%  Program  name  laurant-m 
%  Dieter  R.  Lohrmann  04/26/99 

%  This  program  calculates  azimuth  and  elevation  of  antenna 
%  for  test  on  research  vessel  Lauren  at  NSWC,  Panama  City  FI 
%  scheduled  for  July  15/16  1999 
%  testdata  file  in  testd.m 
%  input  data: 

%  input  is  gps  data  gpsp  of  plane  and  gpss  of  ship, 

%  resulting  in  the  directional  vector  x: 

%  vector  X  points  from  the  ship  to  the  aircraft. 

%  input  data: 

%  gpsplati=  gps  latitude  of  plane  (degrees) 

%  gpsplong  =  gps  longitude  of  plane  (degrees) 

%  Greenwich  meridin  is  zero,  angle 

%  is  counted  negative  in  east  direction 

%  gpsslati  =  gps  latitude  of  ship  (degrees) 

%  gpsslong  =  gps  longitude  of  ship  (degrees) 

%  angle  increases  from  zero  at  G. 

%  towards  east 

%  gpsph  =  gps  height  of  plane  (m) 

%  gpssh  =  gps  height  of  ship  (m) 

%  beara  =  bearing  angle  of  ship  vs  north  (degrees I 
%  pitcha  =pitch  angle  of  ship  (degrees) 

%  rolla  =roll  angle  of  ship  around  its  length  axis, 

%  positive  cw  looking  towards  bow  of  ship  (degrees) 

%  output  data: 

%  azim  =  azimuth  angle  of  the  antenna  beam  in  reference 
%  to  the  boresite  of  the  antenna  array  (degrees) 

%  elev  =  elevation  angle  of  the  antenna  beam  in  reference  to 
%  the  boresite  of  the  array  (degrees) . 

%  intrinsic  data: 

%  antposaz  =  the  azimuth  position  of  the  antenna  on  the  ship 

%  antposel  =  elevation  of  the  boresite  of  the  antenna  array 

%  The  antenna  is  mounted  on  the  ship's  hull;  it  is  not 
%  on  a  stabilized  platform. 

%  The  vector  x  points  from  the  ship  to  the  aircraft. 

%  It  is  calculated  from  the  gps  data  as  follows: 

ax=4e7*cos ( (gpssiati+gpsplati ) *pi/360) * (gpsplong-gpsslong) /360; 
ay= (gpsplati-gpsslati) *4e7/360; 

%  take  into  account  curvature  of  earth: 
az=  (gpsph-gpssh)  -  ( ax'^2+ay^2 )  *pi/  Ael ; 


x= [ax  ay  az] ; 

range=sqrt  (ax'^2+ay^2+az^2 )  ; 

%  transforming  x  to  the  bearing  of  the  ship: 

xp=rlz (x, pi/2-beara*pi/180 ) ; 

%  transform  vector  x  to  pitch  angle: 
x=xp; 

xp=rly (x,-pitcha*pi/180) ; 

%  Tranf ormation  to  roll  angle: 
x=xp; 

xp=rlx(x,rolla*pi/180) ; 

%  transformation  of  x  to  antenna  position: 
x=xp ; 

xp=rlz (x, -antposaz*pi/180 ) ; 

%  transformation  of  x  to  tilt  of  the  antenna: 
x=xp; 

xp=rly (x, -antposel^pi/180) ; 

%  This  is  now  the  vector  of  the  antenna  beam. 

%  The  azimuth  angle  of  the  beam  from  boresite  is  now: 
azim  =  atan2 (xp (2) ,xp (1) ) *180/pi 
%  The  elevation  angle  of  the  beam: 

elev  =  atan2(xp(3) ,sqrt(xp(l)^2+xp{2)^2) )*180/pi 
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%  testd.m 
gpsplati=29; 
gpsplong=-85 • 3 
gpsslati=29 ; 
gpsslong=-85 -  5 
gpsph=150 . ; 
gpssh=5 , ; 
beara=190; 
pitcha=5  ; 
rolla=-5 ; 
antposaz=270; 
antposel=15; 


%  Program  name  rlx.m 

%  Coordinate  transformation  of  threedimensional  kartesian 
%  Coordinate  system.  It  is  rotated  around  x-  axis  ccw  looking 
%  towards  origin  by  angle  alf. 

%  The  unity  vecdtors  of  the  original  system  are  el,e2,e3  . 

%  The  unity  vectors  of  the  transformed  system  are  epl,ep2,ep3. 

%  The  input  row  vector  to  be  transformed  into  the  ep  system  is  x, 

%  the  transformed  row  vector  is  xp  in  the  p  coordinate  system. 

function  xp=rlx (x, alf ) 

el=[l  0  0]; 

e2=[0  1  0]  ; 

e3=[0  0  1]; 

epl=el ; 

ep2=[0  cos (alf)  sin(alf)]; 

ep3=[0  “sin(alf)  cos (alf ) ] ; 

all=sum(epl . -^^el)  ; 
al2=sum(epl .  *e2  )  ; 
al3=sum(epl. *e3) ; 

a21=sum(ep2 . *el) ; 
a22=sum(ep2 . *e2) ; 
a23=sum(ep2 . *e3) ; 

a31=sum(ep3 . *el) ; 
a32=sum(ep3 .  **^62 )  ; 
a33=sum(ep3. *e3) ; 

aa=[all  al2  al3 
a21  a22  a23 
a31  a32  a33] ; 

xp  =  (aa*x ' ) ’ ; 
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%  Program  name  rly.m 

%  Coordinate  transformation  of  threedimensional  kartesian 
%  Coordinate  system.  It  is  rotated  ccw  around  y-  axis  looking 
%  towards  origin^  by  angle  alf. 

%  The  unity  vecdtors  of  the  original  system  are  el,e2^e3  . 

%  The  unity  vectors  of  the  transformed  system  are  epl,ep2,ep3. 

%  The  input  row  vector  to  be  transformed  into  the  ep  system  is  x, 

%  the  transformed  row  vector  is  xp  in  the  p  coordinate  system- 

function  xp=rly (x, alf ) 

el=[l  0  0] ; 

e2=[0  1  0]; 

e3=(0  0  1] ; 

epl=[cos (alf )  0  -sin(alf)]; 

ep2=e2 ; 

ep3= [sin (alf )  0  cos (alf ) ] ; 

all=sum (epl , *el) ; 
al2=sum (epl . *e2 ) ; 
al3=sum(epl. *e3)  ; 

a21=sum(ep2 . *el)  ; 
a22=sum(ep2 . *e2)  ; 
a23=sum(ep2 .  “*^03)  ; 

a31=sum(ep3. *el)  ; 
a32“sum  ( ep3 . *e2 ) ; 
a33=sum(ep3, *e3) ; 

aa=[all  al2  al3 
a21  a22  a23 
a31  a32  a33 ] ; 

xp  =  (aa*x  * )  * ; 
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%  Program  name  rlz,m 

%  Coordinate  transformation  of  threedimensional  kartesian 
%  Coordinate  system.  It  is  rotated  around  z-  axis  ccw  looking 
%  towards  origin,  by  angle  alf. 

%  The  unity  vecdtors  of  the  original  system  are  el,e2,e3  , 

%  The  unity  vectors  of  the  transformed  system  are  epl,ep2,ep3. 

%  The  input  row  vector  to  be  transformed  into  the  ep  system  is  x, 

%  the  transformed  row  vector  is  xp  in  the  p  coordinate  system, 

function  xp=rlz (x, alf ) 

el=[l  0  0] ; 

e2=[0  1  0] ; 

e3=[0  0  1] ; 

epl=[cos (alf )  sin(alf)  0]; 

ep2=[-sin(alf )  cos (alf )  0]  ; 

ep3=e3 ; 

all=sum(epl . *el) ; 
al2=sum(epl, *e2) ; 
al3=sum(epl , *e3) ; 

a21=sum(ep2, *el) ; 
a22=sum(ep2 , *e2) ; 
a23=sum(ep2 . *e3)  ; 

a31=sum(ep3 , *el) ; 
a32=sum(ep3 , *e2 ) ; 
a33=sum(ep3 . *e3)  ; 

aa=[all  al2  al3 
a21  a22  a23 
a31  a32  a33]  ; 

xp  =  (aa*x* )  * ; 
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%  file  name  gztolola.m 

%  calculates  longitude  and  latitude  (degrees) 

%  using  geocentric  coordinates  aax^aay^aaz  (m) 

%  The  apex  of  the  cartesian  right  hand 
%  coordinate  system  is  in  the  center  of  the  earth 
%  The  x-axis  passes  thru  the  Aequator  at  the 
%  greenwich  meridian.  The  Y-axis  is  in  the  aequat. 
%  plane,  and  the  z-axis  goes  thru  the  north  pole. 

%  h  =  heigth  above  surface  of  earth  (m) 
r=sqrt  (aax^2+aay''2+aa2^2)  ; 
la= 18 0/pi* as in { aaz/ r ) 
lo=180/pi*atan2  (aay,aax) 
h=r-4e7/2/pi 
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[1]  Hiitte,  Des  Ingenieurs  Taschenbuch ,  28  th  edition,  p.  1 15  : 


fi  Koordlnatentransformation.  (O;  Cj.  e^,  e^)  und  (O';  c,.  e^)  seien 

zwei  kartesische  Koordinatensysteme.  Zusammenhang  rwiscfaen  ihnen  ist 

feStgelegt  durch  aO  ^  a  ^  ^  . 

ei-<inei  +  <2<ae,+<i<3e3,  ez -aue( -f-aaiC' +a»<ei  (t  =1,2.3). 
a<i -€<  Ck -cos(ei,  Ci). 

(e;.  Ci)  ist  der  Winkel  zwischen  der  x<-Aclise  und  der  Xfc-Achse.  Es  gilt: 
■I  .  fur  1-^1 

P  sei  ein  beUebiger  Punkt  mit  dem  R.V.  t  im  ersten  und  den  R.V.  t'  im 
zweiren  System:  /«/  .  a. 

T  jCjCi  +  XjC,  +  X3C3  ,  V  -  *16,  +  *363  +  XjCs , 

OP  - 1 .  0^  -  CW+OP-u  +  r-t'. 

Die  Koordinaten  von  P  im  ersten  System  hkngen  mit  den  Koordmaten  x't 
von  P  im  zweiten  System  wie  folgt  zusammen: 

xi~a'i+aiixi+<n%x^+aitxz  (t-1.  2.  3). 

Bemerkunff-  Vektor  a  -  O^entspricfat  oiner  Parallel verschiebung  des  Koordinaten* 

.„,JJ:  d"  ^&z  (S  UV  a  -  Jibt  eta.  Dr.l.ua«  d»  Kootd..atan.,s«».  a».  E. 

Ut  Jede  Matrix  mit  dleser  Slgenscbalt  heifit  ortbogonal. 
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